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Abstract 

We use the Cottingham method to calculate the pion and kaon electromag- 
netic mass differences with as few model dependent inputs as possible. The 
constraints of chiral symmetry at low energy, QCD at high energy and exper- 
imental data in between are used in the dispersion relation. We find excellent 
agreement with experiment for the pion mass difference. The kaon mass dif- 
ference exhibits a strong violation of the lowest order prediction of Dashen's 
theorem, in qualitative agreement with several other recent calculations. 



UCTP-9-96 
UMHEP-428 
hep-ph/ 96 11331 



1. Introduction 



The calculation of the electromagnetic mass differences of pions and kaons has 
recently been quite an active topic in the field of chiral perturbation studies. 

This is partially due to the interest in the values of the light quark mass ratios, for 
which we need to be able to separate electromagnetic from quark mass effects 
The realization that Dashen's theorem 0, relating pion and kaon electromagnetic 
mass differences in the limit of vanishing rriu, md, mg, could be significantly violated 
in the real world ^, |10| has added to the importance of a direct calculation of 
these electromagnetic effects. Moreover these calculations have an intrinsic interest 
as state of the art investigations of our ability to handle new types of chiral calcu- 
lations. The classic studies of chiral perturbation theory |ll], |12| are being extended 



to calculations where one must obtain more detailed information of the intermediate 
energy region using dispersion relations (or sometimes models). The electromagnetic 
mass differences are nonleptonic amplitudes which are a challenge to calculate in a 
controlled fashion. It is our goal in this paper to calculate these mass differences as 
well as we can at present. 

Our tool is the Cottingham method for calculating electromagnetic mass differ- 
ences. As explained more fully in Section 3, this converts the mass differences ampli- 
tude into a dispersion integral over the amplitudes for 77r inelastic scattering. As we 
learned in the 1970's from the study of jp inelastic scattering, the physics of such a 
process is reasonably simple. The elastic scattering is well known. At low energies, 
one sees the inelastic production of the low lying resonances. In our study we take 
these resonances and their coupling constants from experimental data. At high ener- 
gies one enters the deep inelastic region for which perturbative QCD can be used. It 
turns out that in the pion mass difference the deep inelastic region cancels out both 
at zeroth and first order in the quark masses. This leaves the mass differences to be 
dominated by the lower energy region. 

There are a series of constraints on the calculation which are important for giving 
us control over our method and results. The most important of these are: 

1. There exists a rigorous result for these mass differences, exact in the limit that 
rriq ^ 0, (g = u, d) which states that in this chiral limit the pion mass difference 
is 



= j dsslns{pv{s) - Pa{s)) 



where py(s) and Pa{.s) are the vector and axial vector spectral functions mea- 
sured in e^e~ annihilation and in r decays []13[. This is a powerful constraint 



because it requires that the full calculation differ from this only by terms of or- 
der or higher, and must reduce to this as — > 0. Many of these deviations 
are kinematic in origin and hence are well tied down by this constraint. 
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Dashen's theorem states that in an SU{3) extensions of this same hmit = 
0,q = u, d, s) that the kaon mass difference is equal to the pion mass difference. 
This means that similar physics enters both amplitudes and one is able to focus 
more directly on SU (3) breaking. 

The low energy structure of the Compton amplitudes 771 — > 771 and 'jK — > 'jK 
are known rigorously from chiral perturbation theory [|T4|-[|l^ and the process 



in the crossed channel 77 —>■ tttt matches well with experiment fl^ 



QCD gives us important information about the high energy behavior of the 
dispersive integral, with the result that Am^ is finite up to order m^, while 
Am|- has at most a logarithmic divergence at order m^, which is to be absorbed 
into the u, d quark masses. This is very useful in pinning down the high energy 
parts of the calculation. 

The medium energy intermediate states are known directly from experiment. 
This region is the most difficult to control purely theoretically, and so we rely 
on experimental data to overcome our inability to provide a first-principles 
theoretical calculation. 



These properties are important ingredients for the reliability of our method. While 
there are still some approximations and educated guesses involved in the matching 
up of the various regions of the calculation, this method is more than just another 
model and represents the real world as well as is possible in analytic calculations at 
present. 

While we estimate that our uncertainty is about 10% for pions, and 20% for kaons, 
our calculated value for the pion mass difference agrees excellently with experiment 
(Amf = 4.54 ± 0.50 MeV vs. Am^^^* = 4.60 MeV). In the case of kaons, our 
calculated value is Am^ = 2.6 ± 0.6 MeV, indicating a strong breaking of Dashen's 
theorem (Am^-^ = l.?)MeV) in agreement with many other recent works 0-0. 

In the next section, we briefly review the physics and history of the calculations 
of electromagnetic mass differences. Section 3 presents the basics of the Cottingham 
method, while Section 4 describes our application of it to the pion mass difference. 
The kaon mass difference is studied in Section 5, and we summarize our findings in 
Section 6. 



2. Review of the Problem 



The mass differences of kaons and pions 



Am^^^'* = m^± - m^o = 4.5936 ± 0.0005 MeV , 

Am^^* = mK± - niKO = -3.995 ± 0.0034 MeV , (2) 

or 
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Ami = 2m^r^'>Am^ = (1.2612 ± 0.0001) 10"^ GeV^ , 

Am^ = 2m^K^^AmK = (-3.9604 ± 0.0035) 10"^ GeV^ , (3) 

where m^K^' = |('^(x,7r)± + "^(x,7r)o), are due to two sources: quark masses and 
electromagnetic interactions. The difference in mass of the up and down quarks 
can produce isospin breaking in hadron masses. However, because the quark mass 
sphtting is AI = 1 and the pion mass difference is only sensitive to AI = 2 effects, 
the pion mass difference only receives contributions of second order, i.e., (m^ — m„)^. 
In fact the leading effect of this order is calculable in chiral perturbation theory 

A 2\ 1 {m^-rndf 2 (.^ 

Amt] = 7fTLt± , 4 

^Vqm A{mu + md){ms-m) ^ ' ^ ^ 

and is quite small. To the level of our approximations we will neglect this quark 
mass effect and treat the pion mass difference as purely electromagnetic. The kaon 
mass difference, on the other hand, does receive an important contribution linear in 

rrid - m„ 

Ami) =^^^^^ml+0{{m^-mdr) ■ (5) 
/ QM rriu + ma ^ ' 

This relation is one of the primary sources of information on quark mass ratios. For 
it to be useful we need to known how much of the kaon mass difference is due to 
electromagnetic interactions. 

We have one handle on the electromagnetic mass differences which comes purely 
from symmetry considerations. The electromagnetic interaction explicitly violates 
chiral SU(3) symmetry, and its effect can be described within the chiral energy ex- 
pansion. At lowest order, which is order the unique effective Lagrangian with the 
right symmetry breaking properties is 

C^ = gEMTr{QVQV^) . (6) 

This Lagrangian produces no shift in the masses of neutral mesons and equal shifts 
for TT"*" and ^ so that it results in 

Ami = ■ (7) 

This equality is known in the literature as Dashen's theorem ||^. It is valid in the 
limit of vanishing quark masses (u, d and s) and hence of massless pions and kaons. 
There are a large number of effective Lagrangians possible with extra derivatives 
and/or factors of the quark masses, so that Dashen's theorem will receive corrections 



of order m^ or equivalently of order mj^ |T8|. Unfortunately the coefficients of the 



higher order Lagrangians are not known, so that one cannot obtain the corrections 
to Dashen's theorem from symmetry considerations. A direct calculation is required. 
In order to obtain the electromagnetic mass shifts, one must calculate 
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= ^ I d'x{n{p) \TJ,{x)M0)\7r{p))D^/{x) (8) 

as in Figure 1. 



Figure 1: Electromagnetic self-energy. 

In momentum space this is 



~ 2 J (27r)4 q^ + ie ' 



T^M\p-q)=^ / d^xe-''^-%'K{p)\TJ^{x)JM\<P)) ■ (10) 



where 



This calculation is different from standard calculations within chiral perturbation 
theory, because we need to be able to explicitly calculate (and not just parametrize) 
the medium energy and high energy contributions. 

There are a few things that we know rigorously about the calculation. Within 
QCD, the high energy renormalization of a quark mass involves a logarithmic di- 
vergence which is proportional to the quark mass itself. Therefore the pion mass 
difference can pick up divergence's only proportional to the second power of the 
quark masses, which will go into defining renormalized masses in Eq. (4). In our 
approximation, or more strictly in the chiral limit, the pion electromagnetic mass 
difference is finite. For the kaon there may appear a divergence of order am„ or am^, 
i.e., suppressed by one power of the light quark masses. This goes into a renormal- 
ization of the quark masses in Eq. (5). In principle there is an ambiguity about 
how much of the electromagnetic interaction goes into the renormalized values of the 
quark masses. This can only be solved by a precise renormalization condition which 
defines the renormalized quark masses. However because this ambiguity is propor- 
tional to airiu and am^, while the kaon mass difference needs only one factor of a or 
{rriii — rriu), this ambiguity is tiny and is far below the sensitivity of our calculation. 

The earliest attempts at explicit calculations (Riazuddin and Socolow pO| ) 



appeared plausible but can now be recognized as mistreating the chiral portions of 



4 



the calculation. The earliest valid method, and still a remarkably beautiful result, 



came in the work of Das et. al. ||2T|. Here soft pion theorems were used to turn the 
matrix element of Eq. (10) into a vacuum polarization function, which in turn can be 
written as a dispersion relation in terms of the spectral functions of the vector and 
axial vector currents, yielding the formula quoted in Eq. (1). Since QCD satisfies the 
chiral and high energy properties assumed in the original derivations, this remains an 
exact statement of QCD in the limit = = 0. The original authors saturated 
the spectral functions by a single vector and axial vector pole satisfying the Weinberg 
sum rules [^, leading to a remarkably good value Am^°* = 5.0 MeV. More recently 
this sum rule has been explored using the measured spectral functions from e^e~ 
annihilation and r decay, plus QCD constraints [0. These show that the physics 
of the pion mass difference is remarkably simple in the chiral limit with the most 
important effects being those of the lightest resonance contributions. The Das et. 
al. calculation remains a benchmark for other calculations and will be an important 
constraint on our work. 

Through the experience of the past decade of studies of chiral perturbation theory, 
we have gained some insight into the physics of intermediate energies. This lead to 
model attempts to calculate electromagnetic mass differences [^. These calculations 
showed a large breaking (up to a factor of 2) of Dashen's theorem due to mass effects. 
To a large extent the violation of Dashen's theorem has a simple kinematic origin in 
the pseudoscalar propagators of the one loop diagram.El Lattice simulations have also 
recently started to be applied to this problem. They also see a significant violation 
of Dashen's theorem, ( Amx = 1-9 MeV ) WU. 



3. The Cottingham method and meson mass shift 

The nonleptonic matrix element which we must calculate is given in Eq. (10). If 
we decompose the Compton amplitude in terms of gauge invariant tensors, we can 
define 



Tf,u{q'^,p- q) = Di^^Ti{q^,p- q) + D2f,uT2{q^,p- q) 

qf^qu 



q' 



D2,u = ^[Pf^-^q,j [Pu-^quj . (11) 

We have used the standard definitions for these tensors. Note that in the soft-pion 
limit, i.e. — > 0, the combination D2fj,uT2 vanishes as we will see in the following 
section. 



^Ref. j|] has an error in one of the mass effects, as described later. We disagree with the 
methodology of a paper which attempted to correct this problem |^ , and agree with the critique of 
101 which is contained in Ml. 
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A first step consists of a rotation in the complex plane and a change of variables. 
We work in the pion rest frame, p-q = rriT^qQ. Since the singularities in T^^ are located 
just below the positive real axis and above the negative real axis in the complex go 
plane, the integration over q^ may be rotated to the imaginary axis, qo — > igo, without 
encountering any singularities. After this transformation the integral involves only 
spacehke moments for photons, i.e., q^ = — = ~{Qo + q'^): the mass shift 
becomes 

A change of variables from {q, go) to (Q^, z^), where v = m-^qo, involves 

d'qdqo ^27C dQ' — ^/m^Q^ _ ^2 ^ (13) 

Jo Jm^Q ^ 

which converts the mass shift to 



5m^ 



dQ^ 



167r^ Jo J-mnQ Q"^ 



roo p7 



+ 



167r^ Jo ' J-runQ 
-3T,{-Q\tiy)+ f 1- 



T2{-Q\iy) 



(14) 



This has reduced the mass shift to an integral over the forward Compton scattering 
amplitude for space-like photons. 

The reduced Compton amplitudes Ti and T2 are presently required to be evalu- 
ated at imaginary momenta, iu. However they can be written in terms of physical 
amplitudes via dispersion relations. The Compton amplitudes are known to obey 
dispersion relations in the v variable with that for Ti requiring one subtraction. 



Ti(q\v) = r,tf,0) + - 



dv'^ ImTi{q\v') 



i,/2 j,/2 _ ^2 



TT Jo y'^ — v"^ 



(15) 



The imaginary part of the forward scattering amplitudes ImTi are defined as 
electron scattering structure functions 



-lm%{-Q^,v) = Wi{-Q^,v) , fori = 1,2 

TT 



(16) 



After employing these dispersion relations, the integral over v can be done explicitly 
with the result 
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+ 30= ^^lV,(-(3^,/)A,(^' 

Jo V \m^Q'^ ^ 



where 



A2(y) = + + + i . (18) 

These manipulations have transformed the mass shifts into integrals over the structure 
functions in the physical region, as well as the subtraction term Ti(— Q^,0). The 
(Q^, z/) plane is shown in Fig. 3.1, as is the physical region where Wi ^ 0. 



4. The pion EM mass difference 



In this section, we describe the details of the calculation of the electromagnetic 
mass difference of the pion. The exact result of the Das et al. |^l]] calculation 
in the chiral limit involves the difference of spectral functions pv{s) — Pa{s)- This 
difference is entirely determined by the leading vector and axial-vector resonances. 
Therefore, our first step is to study the low energy chiral amplitudes supplemented 
by the interactions of vector and axial- vector resonances. These have been previously 
studied in a model field theoretic calculation [Q. We correct a technical mistake in 
that work (which was also noted in [0]), and transform the results into our dispersive 
framework. This allows us to show how the Cottingham method merges with the 
chiral limit result of Das et al. as m.^ — ^ 0. 

We subsequently generalize the calculation by treating the resonances more real- 
istically and adding in other ingredients to the amplitude. The former improvement 
involves the replacement of the "narrow-width" treatment of the resonances, which 
occurs in any field theoretic treatment, by spectral functions which account for the 
energy variation and width of the resonances. To complete the ingredients to the 
calculation, we add resonance transitions not accounted for previously and also the 
deep inelastic continuum. The resonance couplings follow from experiment, and their 
presence in the Compton amplitude is confirmed by the comparison of theory and ex- 
periment in 77 -—>■ vr^Tr*^ |T^, ^ |T^, The deep inelastic region cancels in the mass 
difference to the order that we are working, so that we include only a few comments 
on the matching of low and high energies. 
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Lagrangian with spin-1 resonances 

Our starting point for calculating the pion Compton scattering amplitude is a 
Lagrangian which includes chiral terms and vector {J^'-' = 1 ) and 

axial- vector {J^^ = 1^^) couplings, introduced by Ecker et al. p4|, This La- 
grangian provides an accurate description at low and medium energies (up to ~ 1 
GeV). 

+ ^Tr {V,jn + ^Tr (V.^uV) 

+ ^TT{A,jn , (19) 

The notation is defined in the appendix. The relevant terms after expanding the 
above Lagrangian in terms of pion, photon and spin-1 resonance fields are 



p2 



iGv 



F2 
2eG 



F2 



at 71 



(20) 



The Feynman diagrams which contribute to the pion Compton scattering ampli- 
tude, given by the above Lagrangian, are shown in Fig. ^ It is convenient to classify 
these diagrams in three groups, which correspond with three gauge invariant terms 
of the amplitude. 

The first term encloses the contribution given by Fig. ^.a, and Fig. ^b. 



2Di 



(21) 



2\\2 



+ Ami\GM')\ 



— + 



— {p + qY — it ~ ip ^ l)^ — "^e, 



D 
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(a) 




(e) 



Figure 2: Compton scattering diagrams for the meson resonances, (a) Elastic diagram, 
(b) Pseudoscalar seagull diagram, (c) Vector resonance seagull diagram, (d) Axial-vector 
resonance intermediate state diagram, (e) Pion form factor diagrams. 



where we have used the vector resonance dominance approximation in the pion form 
factor, i.e., 

= (22) 

TT 

This is equivalent to saturating the pion form factor with the rho meson resonance, 

2 

G.(g') = ^^ . (23) 

The second term, which we call the vector seagull, is given in Fig. ^.c. It differs 
from the pseudoscalar seagull because one of the photon lines interacts through a 
vector resonance, 

^2 J2 



TfKq\v^) = -^%^-,B,^.. ■ (24) 

The third group, due to the axial- vector intermediate state, given by the diagrams in 
Fig. ^.d is 

T('.Hg^P■g) = 



\ m\ — (p + qY — ie 

^ {p-q - q'^f + q^ {m\ -jp- qf) \ ^ 

m\ — {p — qY ~" J t^'^ 

F^m\ \m\ — {p + qY — it m\ — {p — qY — it j 



(25) 



The numerical values that we use for the parameters involved in the previous five 
equations are 



= m^± = 0.13956995 ± 0.00000035 GeV , 
= 0.0924 ± 0.0003 GeV , 
mv = 0.7699 ± 0.0008 GeV and 

Fv = 0.1529 ±0.0036 GeV . (26) 
For the numerical results in the soft-pion limit, we use the axial- vector resonance 



parameters, m^, and Fa_, obtained by the Weinberg sum rules p2[ in the narrow- 
width approximation. 
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■,(WSR) 



F2 



F2 = 0.1218 ± 0.0045 GeV 



m 



(WSR) 



Fv my 



0.9664 ± 0.0427 GeV . 



(27) 



The corrections to the final m^, and Fa which we use in our final numerical answer 
at the end of this section, even though necessary to cancel divergences, are minimal 
and also have a minute effect on the numerical results for Am^. 

Beyond the narrow-width approximation 



The above analysis utilizes zero-width p and ai resonances. This "narrow-width" 
approximation is a poor description for these resonances since they are not particu- 
larly narrow, especially the ai. The full resonance spectrum can be taken into account 
by employing the spectral function - Kallen-Lehmann - representation [26, 27 1. Fur- 
thermore, this representation includes the effect of higher mass resonances with the 
same quantum numbers such as p"s in the vector case. The spectral function represen- 
tation of these resonances generalizes the spin-1 resonance propagators encountered 
in the Compton scattering amplitudes, Tj^}{q^,p ■ q), (for i = 1 to 3) given above. 



— — ie 



ds- 



q^ — le 



{2i 



The sum of the three terms, equations (^, |2J, |2^), of the pion forward Compton 
scattering amplitude in the spectral function representation reads 



{2 a 

F| 7o s \ s — {p + qY — ie 



2m2 



+ D2^^ ^Ami\G^{q')\ 



1 



+ (g ^ -q) 



+ 



1 r^^.pA[s) 



ds 



ml — {p + qy — ie 

+ (g ^ -q) 



2 2 

-m^q 



F^ Jo s \s — {p + qy — ie 



(29) 



The narrow-width result can be readily reproduced by letting. 



V,A' ' 



where the values for Fy, Fa, my and are given in Eqs. (pq, 27). 



(30) 



The Cottingham approach 
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Since we have a complete form for the Compton Scattering amphtude we could 
directly calculate the pion EM mass difference with Eq. (|]), avoiding the Cottingham 



method altogether. Nevertheless, the Cottingham method |^ allows us to gain 
control and insight into the calculation. This method requires the break down of the 
scattering amplitude into subtraction and structure functions terms, which are easily 
extracted from Eq. (^), 

2 r, ^ 



V /"OO 

T,{-Q\Q) = -2 + — / dsp^{s 

Ff Jo 



s + Q 



2 



2 Ri ^ Q 




2 




p^Q^ 5{s-p^ + -2u) , (31) 



Before we describe each of these terms in detail, it is useful to be more familiar 
with their domain in the (z/, Q'^) plane. The subtraction term is the value of Ti along 
the negative axis. The structure functions are limited by kinematics to a sector 
of the first quadrant, their domain is better understood if we introduce the Bjorken 
scaling variable, x = The allowed kinematic ranges for the variables involved are: 



0<(5^<C)o;c<z/<oo; and < a; < 1 , where c = — . (32) 

The domain of both structure functions in the (z/, Q^) plane, covers the area in 
the first quadrant which lies between the positive u axis {x = 0) and the elastic line, 
{x = 1). This is the unshaded region shown in Fig. ^for pion kinematics. Within this 
sector, the figure shows other lines of constant-x which help describe the structure 
functions in the scaling region. It also shows two lines of constant s = m\ which mark 
the region where the resonant intermediate states are the dominant contribution. 

The scaling region for nucleons is the region above ~ 1 GeV. It is described 
by perturbative QCD. The relevant degrees of freedom are quarks and gluons, and 
the structure functions are described in terms of quark distribution functions, which 
depend only on x if we neglect logarithmic deviations. In this approximation, the 
structure functions are constant along the constant-x lines. 

The resonance region in the (z/, Q"^) plane is described by the two dashed lines 
parallel to the elastic or a; = 1 line. These lines satisfy the equation for constant 
squared invariant mass of intermediate resonant states, 

Ml = {p + qf =p' + 2v -Q'' . (33) 
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Figure 3: {i^,Q'^) plane. Units in GeV^. Unshaded region is the domain of the structure 
functions. Solid lines are for x = 0, 0.2, 0.4, 0.6, 0.8 and 1.0. Dashed lines are for M|. = 
1.6 and 3.2 GeV\ 



Since the graph is for pion values we choose the ai as the first resonance, nia^ = 1.26 
GeV. We also include a second resonance with mass uir = nia^ in order to show 
the position of possible higher resonances. 

This resonance region is described by chiral Lagrangians which include spin-1 
resonances such as Eq. (|1^), @, |2^. The structure functions obtained through 



the chiral Lagrangians in the narrow-width approximation are constrained to the 
elastic line, the Oi line, plus other parallel lines corresponding to possible higher 
resonances. If we include a finite width for the resonances, these lines become bands 
whose thickness is proportional to the resonance width. 

The usefulness of applying the Cottingham method arises from the breakdown of 
the scattering amplitude into the three terms shown in Eq. (0). We gain control 
because we can make reasonable assumptions and establish constraints on the pion 
structure functions. At the same time, it is possible to relate the subtraction term to 
the soft-pion limit. 

All of the resonances couplings will contain form-factors which suppress the effect 
of an individual resonance as oo. We will assume that the fall-off of all such 

photon form-factors will involve a scale which is a typical vector meson mass. We now 
turn to the procedure to introduce these form-factors in our dispersive framework. 
This has a subtlety in that some naive structures for this form-factor could upset 
the soft-pion limits in our formulas. We will chose a form which is well behaved in 
the soft-pion limit. The form- factor also solves what would appear to be a problem 
in the present inclusion of resonances, i.e. the structure function Wi{Q'^,v) given in 
Eq. ( |3T| ) has terms proportional to v"^ and which would generate divergences for 
large Q^. This is does not occur in the presence of the form-factors. This divergent 
behavior is clear if we calculate v) along the lines of constant s = mj^ . We 

use the 
obtain 



use the delta function to eliminate the u dependence, i.e., u = |(s — p'^ + Q"^) , to 



Wi-''\Q\ V) = Us-p'- Qr , (34) 
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which diverges as Q^. A structure function for a given resonant state cannot diverge 
for large without violating unitarity. For a single resonance, as it is the case 
in question, the structure function must go to zero if we follow a line of constant 
invariant mass, s, to high energies. In order to achieve this behavior, we introduce a 
multiplicative factor which forces its convergence. This factor, K{v^ s) resembles the 
form factor obtained for the elastic term through the vector meson dominance model. 
In addition, it has the following properties. 



lim K(z/,s)(^^*-^) ~ — , 



Q2=0 



1 . 



(35) 



These properties ensure that the subtraction term is left unchanged, and that the 
structure function will converge for large Q^. The form factor is normalized in order 
to agree with the previous result at = for fixed s. The form factor that satisfies 
these conditions is 



my 



m 



1v 

1 + r7— 

s 



(36) 



where 



s — p'^ 



1 + 



P 



1 



(37) 



and = . We have chosen the appropriate value for the vector meson mass, 
my = nip . The factor K{i', s) above, also has the property of being very close to the 
p contribution to the pion EM form factor for s = p^, as it can be seen in Fig. |. 

The inclusion of this factor in our analysis is easily achieved through the substi- 
tution 



The structure functions and subtraction terms read 



(3^ 



Tl(-g^o) 



Ri 



F^Jo 



dsp 



Ri 



5 _ p2 _|_ Q2 



S 



-Ti JO s 
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1 2 3 4 5 



Figure 4: Factor /^(z/, s) and pion EM form factor. Pion EM form factor, solid line. 
K{u,ml), dashed line. K{u,l GeV^), dash-dotted line. K{u,2 GeV^), dotted line. 
Horizontal scale represents in GeV^. 

^1 r^^^pM. K{u, s) p'Q^ 6{s -p' + Q'- 2v) . (39) 

JQ s 

The pion EM mass difference for the above functions is readily obtained with 
Eq. (0). We choose to break it into terms corresponding to those shown in the above 
equation with an extra subdivision of the W2 contribution which isolates the elastic 
term 



Aml.{Subtr.) = 

AmliElast.) = 
Aml{W2) = 




where = s — + Q^, and Aj(?/) (for i = 1,2 ) are defined in Eq. (pISl). 

High energy constraints 

Even though, we can obtain an explicit formula for the pion EM mass difference 
by adding all the contributions in Eqs. (|iO|), it is necessary to analyze the upper limit 
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of the integral. Adding all the different contributions in Eqs. P0|), and expanding 
the integrand in powers of we obtain, 



+ 



/ 3a \ 



fOO 

-Fl + ds 
Jo 



ds s 



p^(^)-p5(^)(i- 




where, 



mz 



If the first two terms are not zero, they originate linear and loga- 



rithmic divergences respectively. In order to obtain a finite pion EM mass difference 
we cancel them explicitly generating two constraint equations. 



ds s 



ds 



P%s) 1 



2 4 ^ 

s s^ 



(42) 
(43) 



These high constraints become the Weinberg sum rules in the soft pion limit, 
which in the above equations is obtained by letting = 0. We will see later a more 
detailed explanation of this limit and its relation to the subtraction term. 

Due to the introduction of the convergence factor K{i>,s), the above divergences 
originate only from the subtraction term. We incorporate the high constraints in 
the subtraction term of the pion EM mass difference Eq. ( ^OD in order to make all 
the contributions finite. The following procedure removes both divergences. 

Subtract the linear divergence from the subtraction term by means of Eq. (^2]), 



Am1{Subtr.) 



oo Q2 

"^'p^^'^VTq-^ 



+ 



5 _ p2 _|_ Q2 



POO 

/ dsp 
Jo 

/•oo roc I irp' 

F'^ - I dsp^is) + dsp^is) i^l-^ 




(44) 



Integrate over Q"^ and cancel the correspondent logarithmic divergence by subtracting 
Eq. (ID multiplied by In Aq2 , 
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Am'l{Subtr.] 



32^00 47r Fi Uo 



A02 + s 



JO \ 5 / S — p"^ 

Add the terms and take the hmit Aq2 00, to obtain 



(45) 



Am^^Suhtr.) 



a 3 



47rF2 



ds 



s\ns p^{s)-{s-p')\n{s-p')p%s)\l 



2\ 



J9^ 



(46) 



The above contribution is free of divergences. Furthermore, in the soft-pion hmit, i.e. 



= 0, it is equivalent to the Das et al. calculation Finally, we have a useful 
formula to calculate Sm^J^^^^ free of divergences, which was mainly the product of the 
Lagrangian introducing the chiral couplings of the spin-1 resonances. We proceed to 
show the close relation of the subtraction term and the soft-pion limit and to see how 
we reproduce the results of O(p^) chiral perturbation theory with the above scattering 
amplitude. 



Soft-pion limit and its relation to the subtraction term 

In the following we will show that the subtraction term is given by the soft- 
pion limit up to corrections of order p^. In this discussion we refer to non-contact 
contributions as all contributions except the pion seagull term. This term is the only 
one which has both photons interacting at the same vertex and therefore we treat it 
differently in the following discussion. 

The non-contact contributions to the Compton scattering amplitude have the form 



T^:!''\q',p-q) = I jdxe-^^-{T,\T(y,{x)VM) • (47) 
Consider the soft-pion theorem [p9[] . 



hm (7r'=(p)/5|0|a) = [Ql 0\ \a) , (48) 

where (3 and a are arbitrary states and = Jd^xA^^x) is an axial charge. We also 
need the commutators 
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The result of applying the soft-pion theorem to Eq. (|^) is 



(49) 



,-2 



lim T(^^)(g^p■g) = -i d\e^^^ — {^\T {V,^{x)V^ {Q) - A^{x)Aim |0) . (50) 



The two-current time ordered products are related to the spectral functions by 



(0|T(K''(x)n'^(0))|0) = i6atl dspv{s){ng^,-d^dnj 

(27r)4 P + ie 
d^k e-'^'' 
(27r)4 F - s + ze 



d'^k 



-ikx 



(27r)4 F - s + ' 



(o|t(A(^(x)a;;(o))|o) 



(51) 



Upon combining Eqs. (|50D and (0), integrating over d'^x, and using the resulting 6 
function to integrate over d'^k, we obtain 



-2 



g2 + ze 



(52) 



(py(s) -Pa(s)) 



g2 — s + ze 



The contact term or pion seagull contribution to the pion Compton scattering ampli- 
tude is 

T^'^\q\p-q) = 2g^u , (53) 

which remains unchanged in the soft-pion limit. Adding both contributions to the 
pion Compton scattering amplitude in the soft-pion limit, we obtain 



{2 f°° 
-2 + — / ds{p 
Jo 



v{s) - Pa{s)) 



— s 



(54) 



An alternative way of reproducing the soft-pion limit result above, is letting 
in Eq. (^). In order to implement this limit the following relations are useful. 



lim D^ 

lim rP'Do 

lim Ti(o^, p-q) 
Pm^o '-^ ''^ 



, 

Tl(g^O) 



p2=0 



(55) 
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From these relations it follows that the only surviving term in this limit is the sub- 
traction term Ti(g^, 0), 



\im T {q^,p-q) = Di lim Ti(g^O) 



This gives the same result as Eq. (Q) when we identify 

hm (p^is) - p%s)) = (pvis) - pAis)) . 
We can now calculate the soft-pion limit to the pion EM mass difference, 



lim Am'i 



Q2 

ds {pv{s) - pa{s)) 
+ s 



(56) 



(57) 



(5^ 



where = —q^. We follow the procedure described earlier in order to cancel the 
linear and logarithmic divergences occurring in the above equation. The cancellation 
of these divergences imposed by the finiteness of the pion EM mass difference requires 



/■oo 

/ ds{pv{s) - pa{s)) 
Jo 

dss {pv{s) - Pa{s)) 



(59) 
(60) 



These are Weinberg sum rules p2[, obtained in our case as a consequence of the 
finiteness of 6m^J^^'^^ in the soft-pion limit. Subtracting the linear and logarithmic 
divergences in the same way as for the subtraction term, Eqs. (|4^ - p6D , we obtain 



... 2 ci; 3 
lim Am^ = - — — 

Pm^o Arc F2 JO 



^oo 

/ dss\ns{pv{s) - pa{s)) 
Jo 



(61) 



This is the result obtained by Das et al. 



Finally, we evaluate Eq. (|6l| ) in the narrow-width approximation to obtain the 
numerical result 



lim Ami^^^ = 4.685 MeV 

Pm^O 



(62) 



Spectral functions 

We seek an improved description of the physics of the resonance region with 
the spectral functions Py(s) and p^(s) replacing the narrow- width description. The 
ingredients to the spectral functions clearly are the same resonance states that are 
revealed by the usual vector and axial- vector spectral functions pv{s) and Pa{s). In 
addition we have just seen that in the soft-pion limit there is an exact correspondence 
Py('S) — p5(s) = pv{s) — Pa{s)- This leads us to utilize the experimental spectral 
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functions determined in reference |T^ in order to produce a shape for py{s) and 
p^(s). In both the soft-pion hmit and the full Cottingham calculation at order m^, 
the high energy continuum cancels in the mass shift. We therefore separate each 
spectral function into two contributions, one due to the resonances and the other due 
to the high energy continuum common to both vector and axial-vector channels. The 
resonant part is chosen to match the resonances revealed in the phenomenological 
analysis of the data in [T^. These spectral functions are then slightly altered to obey 
the full constraint equations including terms of Eqs. ( pSj) and ([86|). A continuum 
contribution, common to both vector and axial- vector channels, was included in ||13| , 
but is here kept separate from the resonances. The result of this is that we identify 



PvA'V 

with a continuum contribution 



(63) 



P^Ais) 



(64) 



The precise identification of the continuum is not unique, but since the difference 
of spectral functions enters, reasonable variations do not produce a large final effect. 
The specific form that we use is shown in Fig. 5. 



0.06 



0.04 



0.02 




Figure 5: Vector and axial-vector spectral functions. The graph shows (V) Py(s), (A) 
p^(s), and (C) pc(s) versus s. The s scale is given in GeV^. 



It is clear that the greatest source of model dependence in our calculation comes 
form the numerical identification described above. Our procedure in setting up the 
calculation in the Cottingham method is very general. However, we do not have 
directly available the experimental structure functions for photons scattering off of 
pions. We have used an identification which is valid in the chiral limit in order 
to provide this numerical input. There could be shifts in the couplings of these 
resonances which are of order m^. These could provide changes in the final answer 
at order which would be of interest to us. This is partially relieved by the fact 



that the analysis of [|T3| was carried out with real world data, not strictly in the chiral 
limit. Thus the masses, widths and shapes of the resonances will accurately reflect 
physics with ^ 0. Likewise we know that in the narrow- width approximation we 
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have the right description, so that we don't see a source of major uncertainty due to 
the non-zero widths. This means that our model dependence comes from possible 
dependences in the resonance couplings, and our implicit assumption is that these 
are smaller than the dependence from the propagators. We have not been able to 
find a way to do better than this in the phenomenological analysis. 

Comparison 0{E^) Chiral Perturbation Theory 

We also like to compare our method to the standard chiral perturbation approach. 
The lowest energy region of the pion structure function can be described by the chiral 



SU(3) Lagrangian to order p , originally developed by Gasser and Leutwyler 0, |T2 
Besides the elastic and seagull terms, and ignoring pion loops, the only relevant terms 
involved in the pion Compton scattering amplitude are the Lg and Lio terms. 



+ LioTr (u^F^'UFl^,) + other . (65) 



The pion forward Compton scattering amplitude resulting from this Lagrangian 
was calculated by Bijnens and Cornet [0], and Donoghue et al. [|l5l. Their result. 



up to pion loop contributions which are small, is 



-2^1.^ + + loops , (66) 

Expanding our narrow- width result in powers of external momenta and g^, we 
obtain, 

■if ^ 2. n -mlq^ \ (67) 




+ higher order terms in {q'^,p-q,p'^) ■ 
The relations for the Lg and Liq in terms of the spin-1 resonance parameters are 



obtained by inspection from Eqs. (^) and (pTf) . 
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This result is in agreement with Ecker et ah The above equation for Liq is also 
the narrow-width approximation for the sum rule (WO) in reference 
Substituting the narrow-width parameters in Eq. (|^) we obtain, 



Lg = (7.20 ±0.05) X 10"^ , 

Lio = -(5.89 ± 0.65) X 10"^ . (69) 
These are seen to be within reasonable agreement with the experimental values, 

Lg = (7.1 ±0.3) X 10"^ , 

Lio = -(6.84 ± 0.3) X 10"^ . (70) 

The difference between the values in Eq. (^9]) and (|70D gives an estimate for the loop 
contributions which we neglected in the 0{p'^) chiral Lagrangian calculation of the 
Compton scattering amplitude. Besides the loop corrections, the difference can also 
be due to the inaccuracy of the narrow-width approximation. 



Scaling region 

The low and intermediate energy regions of the structure functions are described 
above. To complete the analysis of the structure functions we need to describe the 
scaling region at large values of (z^, Q^). The ingredients and general behavior in 
this region are well known. The structure functions become largely functions of the 
Bjorken scaling variable x = (5^/2z/, with logarithmic variations predictable by 
QCD pO[. This is easy to build into the Cottingham analysis [0. However there is 



not a need to describe the details here since the scaling region cancels in the difference 
between the charged and neutral pions masses, to the order that we are working here. 

In the limit that the u and d quark masses are equal, the deep inelastic structure 
functions of the neutral and charged pions are equal. This leads to 

AitIt, (Scaling) = . (71) 

To the extent that the u, d masses are different, the structure functions may differ. 
However we are calculating the electromagnetic effect in the limit = so that 
we are not sensitive to this effect. 
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V TT'-f contribution 

To have a more complete phenomenological description of the Compton scattering 
amphtude we also include the effect of intermediate vector meson diagrams shown in 
Fig. p. The motivation for introducing these diagrams is the experimental observation 
of the radiative meson decays uj and p —>■ iT'-f, and (p —>■ Tf-f. The effective 

Lagrangian which includes the Yttj vertices is 

£ = e^e^'^^^F^.T^^^TT , (72) 

This Lagrangian is invariant under parity and charge conjugation transformations, as 
well as under chiral rotations. The choice of including the EM field strength tensor 
ensures gauge invariance, and the pion momentum dependence corresponds to the 
correct soft-pion limit for the vertex. 



7 7 7 7 




Figure 6: Intermediate vector diagrams. 

We introduce the spectral functions gvis) to describe the intermediate states in 
Fig. 1^. The normalization of these functions is chosen in order to make the subtraction 
term contribution compatible with the ones obtained for the axial-vector case. The 
narrow- width approximation for gv{s) is 

gvis) = H^6{s-ml) = F^RvS{s - ml) . (73) 

The subtraction term and the structure functions for the intermediate vector meson 
diagrams follow from the Lagrangian in Eq. (|T2|), 

1 

Wi{-Q^,u) = — dsgv{s)K{iy,s){iy^+p^Q^)6{s-p^ + Q^-2u) , 
1 

W^{-Q\v) = — dsgv{s)K{u,s)p'QH{s-p' + -2u) , (74) 
Ft Jo 
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where s) is the factor defined in Eq. (|36|). The factor K{h',s) ensures the 
convergence of the structure functions as in the intermediate axial-vector state case. 
We only need to find the spectral function gv{s) in order to determine the above 
functions. 

There are four possible vector intermediate states for the pion Compton amplitude, 
the p"^ for the charged pions, and the p°, uj and (p for the neutral pion. The coupling 
constants, Rp±, RpO, R^, and -R^ can be extracted from the radiative decays of these 
vector mesons. We refer the reader to references [ITTl E3l B2l RSl for a review and 



examples of obtaining such couplings. The couplings Rp± and RpO are the same if we 
take isospin to be an exact symmetry. This means that the charged and neutral pion 
EM self energies due to the p intermediate state would cancel in the pion EM mass 
difference. However, the u and intermediate state contributions do not present 
such a cancellation. Since isospin breaking effects are generally of small magnitude, 
we shall neglect the intermediate p contribution to the pion EM mass difference. 

We can determine the u coupling, R^^, from the experimental measurement of the 
radiative decay u 7r°7, 

24 nv^ 

Ru. = — , ^ "p., r^^^^ = 5.40 ± 0.32 GeV-2 . (75) 
a [m^ — m^)-^ 

Likewise we determine the coupling, i?^, 

= — , 2 ^ r^^^-, = 0.019 ± 0.002 GeV-2 , (76) 
a (m^ — mj^y 

where we have used the experimental values listed by the Particle Data Group p^ . 
We do not consider the vector meson intermediate state further because its coupling 
is an order of magnitude smaller than the experimental uncertainty of the u coupling. 

We are now ready to determine the spectral function gv{s)- Since the only reso- 
nance involved is the u, we can safely use the narrow- width approximation of Eq. (|73|). 
The width of the u is only 1% of its mass. This is in contrast with the p and Oi res- 
onances for which the widths are 20% and 33% of their mass respectively. In the 
narrow- width approximation we only need m^, taken from p^], and Hy, given by 



Hv = F^^Rv = 0.215 ±0.013 . (77) 

We should be careful when comparing gv{s) to Pv,a{s) since they have different 
units. The relationship among these structure functions will become clear in the 
following subsection. 

The intermediate vector meson subtraction term and structure function contri- 
butions to the pion EM mass difference are obtained by combining Eq.s (|T^ and 
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X 



VQ2 



where A"^ = s — p"^ + Q'^, p"^ = m^, and the functions Aj(?/), for i = 1,2 are defined 



in Eq. (|18]). The extra minus sign appears because the vector intermediate state 
diagrams contribute to the neutral pion EM self energy. 



General treatment of other possible contributions 

At this point we have a fairly complete calculation of the pion EM mass difference 
broken down into different contributions. We have included the spin-1 resonances 
through their lowest order chiral couplings, the scaling, and the intermediate vector 
resonance contributions to the pion EM mass difference. By analogy with the nucleon 
structure functions, we are comfortable with our estimates of the structure function 
contributions. These are small, and even a factor of two correction would amount 
to a small correction to the total mass difference. Therefore, we concentrate in the 
subtraction term contribution estimate. 

The subtraction term has been obtained by calculating the pion Compton scat- 
tering amplitude with the effective chiral Lagrangian for the vector and axial-vector 
resonances of Eq. (|1^) and the effective Lagrangian for the intermediate vector meson 
contribution of Eq. (0). In general, there could be other possible contributions to the 
subtraction term. These could be introduced by higher order effective Lagrangians. 
Their contributions to 6m^^ would be small since they would be of higher order in 
the external momenta p"^ and g^. The terms with higher powers of p"^ are naturally 
small, otherwise terms of higher order in are in principle divergent. The finiteness 
of 6m^^'^ requires that all the higher powers of cancel in the same way that the 
order 1 and cancel due to the Weinberg sum rules in the soft-pion limit case. 

We include all other possible contributions to the subtraction term, not yet ac- 
counted for in the previous analysis, by introducing the remainder term, 

— / dsR{Q',p',s) . (79) 
Jo 

The purpose of including this term is to show explicitly the effect of possible correc- 
tions to our current scattering amplitude and its role in the high energy constraints 
and final formula for the EM mass difference. 

There are some conditions required upon this remainder term. It cannot alter our 
previous soft-pion limit result, therefore, 

POO 

lim / dsR{Q'^,p'^,s) = . (80) 
pm^o Jo 

It is also convenient to use the following notation for its expansion in powers of 
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We have explicitly introduced a factor of in order to make sure that this term 
vanishes in the soft-pion limit as given in Eq. (^). This limit also requires that the 
functions hi{p'^,s) (for i = 1,2) do not have a pole at p"^ = 0. The above equation 
is not a formal expansion in orders of since we have introduced the function 

f{p^, s) in the denominator of the second term. This has been done in order to make 
the integral of this term convergent at low Q^. The reason for choosing the above 
notation will be clear in the following extraction of the subtraction term contribution 
to (5mf^^ 

We rewrite the subtraction term contribution, 

AmliSubtr.) = ^ dQ^T,{-Q\Q) . (82) 
87r Jo 

The subtraction term including the remainder part, except its 0{1/Q'^) contributions 
not present in h^, is 



2 ( 



-PAis) 1^ fl - -) - 9vis)^^ (83) 

s — p^ + Q'^ \ s J s — p"^ + Q"^ 

, 2, / 2 X , P^h2{p^,s) 

+p hi{p , s) + 



We expand the above equation in powers of l/Q^ in order to obtain 



Tl(-g^o) 



2 r'^ 



poo 

+ / ds 
Jo 



sp^is) - ^-^-^Pl{s)-p\s-p')gy{s)+p'h,{p\s) 



1 



+0 



The finiteness of Sm^ requires the cancellation of the linear and logarithmic diver- 
gences, resulting in the constraints. 
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ds 



Pv{s)-p%s) 



p 



^5) 



ds 



2p5(s) - sgv{s] 



9v{s) 



+p^/l2(p^ s) 



. 



^6) 



These constraints also reduce to the Weinberg sum rules when we let = 0. The role 
of the functions hi{p'^, s) (for i = 1, 2) is to include all other possible contributions. The 
above constraints must be satisfied exactly, otherwise the pion EM mass difference 
would be divergent. 

We can use the Weinberg sum rules, Eqs. (|59| ) and (|60|) , to further simplify the 
previous equations. 



POO 

/ ds 

Jo 



gvis) + /ll(p^s) 



p'J^ dsU2p'J,{s)-sgvis] 



+ p" 



9v{s) - 



pii^) 



+ h2{p\s) 



, (87) 
. (88) 



We can use solutions available for the functions p2(s) |jT3[ and gvis) to estimate the 
integrals for the remainder terms. 



roc roc 

/ dsp^hi{p^,s) = / dsp^ 
Jo Jo 



9v[s) 



Pa^s) 



3.0 X 10" 



dsp h2{p , s 







ds\p^ 2p^{s) - sgvis) +p 



9v{s) - 



Pa{s) 



6.1 X 10 



-4 



^9) 



(90) 



As expected, these values are small when compared to the integrals involving Py{s) 
which are the larger terms in the constraint equations. 



roc 

/ dsp^(s) = 3.83 X 10 
Jo 

roc 

/ dssp?;{s) = 5.62 x 10 
Jo 



(91) 
(92) 



The remainder term R{Q^, p"^, s) allows us to satisfy the constraints exactly since it 
introduces a small correction to the previous constraint equations. We can proceed to 
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find the subtraction term contribution by following the steps that we used previously 
in order to obtain Eqs. (^4|)-(^) 



Aml{Subtr.) = 



Even though the functions f{p^,s) and h2{p^,s) are undetermined, we have seen in 
equation (BO) that their contributions to the constraint Eq. (RUI) are small. 



Numerical result 



The total pion electromagnetic mass difference is given by the addition of the 
elastic term of Eqs. ( pOD , the structure function terms of Eqs. (^) and (|78|), and the 
subtraction constant term of Eq. (|93D . The results for the narrow-width approxima- 
tion and for the corresponding spectral functions is given in Table |1]. 



Table 1: Am^ results. 



Narrow-width p^{s), py{s) 





(MeV) 


(MeV) 


Subtr. 


4.306 


4.124 


Elastic 


0.500 


0.500 


Str. Fn. ai int. st. 


0.028 


0.041 


Str. Fn. uj int. st. 


-0.127 


-0.127 


Total calculated 


4.707 


4.538 


Experiment 


4.594 


4.594 



We see from the results that the dominant contribution comes from the subtraction 
term, which is largely the effect of vector and axial-vector resonances, with modest 
dependence on = m^. The elastic term gives the only other significant contribution. 
The modification due to non-zero width is also not large. The overall result is in 
excellent agreement with experiment. 

5. The kaon EM mass difference 



Having set up and tested our methodology for the pion, we now proceed to the 
calculation of the kaon electromagnetic mass difference. The most important effect is 
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that the larger mass of the kaon leads to kinematic corrections in the various formulae. 
There are also changes in the mass, width and couplings of the resonances which we 
extract from the data 

The kaon calculation is very similar to the pion one described in the previous 
section, therefore we will concentrate on the differences that arise in the kaon case. 
The kaon counterpart for the Lagrangian of Eq. (^0]) expanded in terms of kaon, 
photon, and spin-1 resonance fields is 



L = ieA^iK+d^K- - K-d^K+) +e^A''Af,K+K- 
e 



+ 



iG 



K 



ieFA 



2F, 



K 



(94) 



where we have used ideal mixing for the vector meson resonances. 

The major difference between the kaon and pion Lagrangians - Eqs. (20) and 
respectively - is that in the kaon case all the three nonet vector resonances contribute. 
Another difference is that there is an elastic contribution to the neutral kaon self- 
energy. This contribution vanishes in the ideal mixing approximation together with 
the limit where all the vector resonance masses are equal. 

The contribution to the kaon Compton scattering amplitude given by the Feynman 
diagrams of Figs. Q.a and |3.b is 



X 



{GK^{q')f-{GK^{q')f 



m\; — (p + qY — ie 



+ {q^ -q) 



(95) 



where 



and 



du 7TSK+fi(u) , (96) 

u — q^ 



6kMu) = ±lsiu-ml) + ^6{u~ml) + ^6{u-ml) . (97) 
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We have subtracted the neutral kaon contribution in order to be able to use this 
equation in the following kaon EM mass difference formulas. 
The vector seagull contribution, Fig. 0.c, is 



' [du — — 

'rJ w - r 

Finally, the axial-vector resonance intermediate state contribution. Fig. ^.d, is 



+ {q^-q) \d,^ 



Flml^^y m\^^-{p + qy-ie ^' j '''' (99) 

Fr^k^^ - (p + g)2 - ^ ^ ^ J ^'"^ ' 

The axial vector intermediate state for the kaon is less straightforward than for the 
pion since the axial- vector meson Ki^ is an ill-determined mixture of the physical 
states ii"i(1270) and ii"i(1400) [0. We will treat this issue later when we estimate 
the spectral function for the axial-vector intermediate state. 

The breakdown into structure functions and subtraction term of the Compton 
scattering amplitude given by the above three terms T'^^^^\ (for i = 1 to 3), in the 
spectral function representation is 



2 



s + Q 



2 



2 /-^ Ki^^^^_(, P'^ 



Jo ' ' S - p2 + Q2 \^ 3 y ' 

1 f.P^i^) 



Wi{-Q\u) = fds^^{u-Qyk{u,s)6{s-p^ + Q^-2u) 

r J S 



K 

W2{-Q\v) = 4m%6{Q'-2u)[{GKA-Q")y~{GKo{~Q'')f 

^2 

K • 



Ids P^Q^k{iy, s)6{s -p^ + Q^- 2v) , (100) 

S J S 



where u = p-q = rnxqa, and p"^ = Tn\. We have also introduced the convergence 
factor 



K{v,s) = jdu[^-^y {l+f,{u)—]5K+{u) , (101) 
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where 



s — p"^ 



1 + 



s — p 



u 



(102) 



and p^ 



m 



K 



The definition of the kaon vector spectral function is 



Pvis) 



(103) 



where pf{s) is the spectral function introduced in Fig. ^, and the other two are due 
to the and u intermediate states. For these last two it is appropriate to use the 
narrow-width approximation. 

We also include the VP7 vertices in the same way as we did for the pions. The 
effective Lagrangian that we use for the K*K'y vertex is 



K* 



(104) 



The subtraction term and structure functions that follow from the above Lagrangian, 
including the K'^ functions with an extra minus sign to be able to insert them directly 
in the kaon EM mass difference formula are 



Tl(-Q^O) 



2 

r !<- Jo 



p'Q' 



S - p2 + Q2 



dsg^ (s) {u' + p'Q') K{iy, 3)6(3 - p' + Q' - 2p) 



/ dsg^ is)p'Q' K{u, 3)5(3 - p' + Q' 



2u) 



(105) 



The narrow-width approximation is justified in this case due to the small width of 
the K* intermediate states. Therefore, we use the following definition. 



9vis) 



E H'j^,, 6(3 - mj^,o) - Hf^,+ 6(3 - mj^^+) . (106) 

Unlike in the pion case, there is an intermediate vector meson contribution for the 
neutral kaon as well as for the charged kaon, in the SU(3) limit. For the pions, the 
intermediate p meson contribution canceled in the SU(2) limit, (since charged and 
neutral p couplings become equal), leaving only the intermediate u and (p contribu- 
tions to the neutral pion Compton scattering amplitude. 

We determine the K* couplings from the radiative decays K* — > Kj , 



R 



K* 



24 ^ 

a (m^,+ — m|-)3 ^* ^^'^1 



0.70 ±0.06 GeV" 



24 m-" 

1 1^*0 



a (mj^.o - mj^) 



-K-^^KO-y = 1.61 ±0.14 GeV-' , 
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from which we obtain the values for 



Hj^,+ = Fk^Rk*+ = 0.093 ±0.004 , 

Hj^,o = Fk^R^o = 0.143 ±0.006 . (107) 

In the kaon case, the mass difference need not be finite because there can be 
divergences which are absorbed into the renormahzed masses of the up and down 
quarks. However this effect is relatively small because it is proportional to arriu or 
arriii compared to the dominant electromagnetic mass-shift which is simply of order 
a. We assume that the renormalization of the up and down quark masses has been 
carried out, although the precise renormalization prescription is hard to define because 
of the small size of this effect. The remaining electromagnetic effects are finite. 

We can now determine the full high constraints for a finite kaon EM self 
energy given by combining Eqs. ( |100| ) and ( |105| ), and including the remainder terms 
introduced in Eq. (|8T|) , 



ds 



Ki 



±p" 



9'^{s)+hnp\s) 



K (J2 



Fl , (108) 



ds 



9vis) 



, (109) 



where = mj^. These equations are similar to the ones encountered in the pion case. 

Using the above constraints and all available data for the vector spectral function, 
Pv{s) and the narrow-width approximation for gv{s), and the uj, (p, ^^(1420), u;(1600), 
and 0(1680) resonances we obtain for hi = h2 = 



Jdspfis) = 0.0170(12) GeV^ , (110) 

/ 2\ 2 

Jdssp^is) il-^j = 0.0406(16) GeV^ , (111) 

where = mj^. 

We try first the narrow-width approximation to the axial- vector spectral function 



Pa is) 



A (NW) 



cos^ 6k5{s — m]^^ 



(1400) > 



± sin^ 9k5{s 



m 



i<'i(1270). 



± Rl5{s 



m 



(112) 
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This parametrization includes the i^^i(1270) and 7^1(1400) resonances as well as a 
higher mass resonance Ki{HR). The input parameters to for this calculation are 



TT 

yR = — and mKi{HR) = 2.0 GeV . (113) 
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These choices are sensible but arbitrary. They fix the values of and Rk through 
the constraint Eqs. (|110|) and ( |111D . The obtained values for and Rk show a 



sizeable dependence on the choice of itik-^^i^hr)- The results for the different contribu- 
tions to 5m^*^ obtained by this narrow-width approximation are given in Table 
In particular we find that the subtraction term contribution is very large. How- 
ever, this contribution varies from 5m^"''*^'' ~ 2.3 MeV for rriKj^iHR) = 1-8 GeV to 
5mg"'*"^ -3.1 MeV for "rriKiiHR) = 2.4 GeV. These numerical results only constitute 
a very rough estimate. This is already indicated by the large dependence in ^/^^(///j), 
and once again, it involves the narrow-width approximation for broad resonances. 
There is another constraint on the axial-vector spectral function due to r decay 




where 



|2 

= ' = 0.6633 GeV-2 . (115) 

The data for the tau lepton decay, r UrKTnr , gives the branching ratios p5l , 



BR{t- ^ UrK-7i+7i-) = (0.40 ±0.09)% , (116) 

BR{r- ^ UrK\'n^) = (0.41 ± 0.07)% , (117) 
BR{t- ^ i^^K-n\^) = (0.09 ±0.03)% , (118) 
BR{t ^ UrKnn) = (0.90 ±0.12)% . (119) 

The last branching ratio is the sum of the three different decay channels with the 
uncertainties added in quadrature. Even though we expect these branching ratios to 
be dominated by the axial- vector channels specially the i^'i(1270) and the A'i(1400), 
there should also be a contribution due to the vector resonance i^'*(1410) This 
resonance will contribute through the decay process i^'*(1410) — >■ i^*(892)7r Ktttt. 
The branching ratio for i^'*(1410) /C* (892)77 is greater than 40% at 95% confidence 
level 0, and BR{K* (892) Ktx) ~ 100% . Therefore, the tau branching ratio into 



the strange axial- vector channels should be somewhat lower than stated in Eq. ( |119| ). 

We obtain the shape of (s), up to rriK-Kn = 2.1 GeV, from diffractive production 
experimental data obtained by the ACCMOR collaboration in 1981 [^. The data was 
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extracted from 200,000 examples of the reaction K^p — * K^n^Ti^p. The intensity 
for the 1+ channel gives us the shape of p^(s). 

We add a high energy tail to the data, up to rriKn-K = 2.54 GeV, which decreases 
quadratically. The constraint equations favor this quadratic choice instead of other 
simple parametrizations. The final normalization of the spectral function is obtained 
by enforcing the constraint Eqs. ( |110| ) and (|1 1 1| ) . If we define 



A Pa[s) , where j dsp^is) = 1 
we obtain = 0.144 GeV, and the p^(s) given in Fig. 0. 



(120) 




Figure 7: p^'^\s) vs. mKT,j, (GeV), for i = 1 to 3 . 

The choice of results in an estimate for the tau branching ratio for the i^*(1410) 
vectorial resonance, 

BR{t ^ UrK*{U10)) = (0.46 ±0.13)% . (121) 

this value could be extracted from data by doing an angular momentum analysis of 
the Ktitt final state. 

The final results for the different contributions to Aitlk are given in Table |[ We 
estimate the uncertainty of the total kaon EM mass difference to be a{AmK) ~ 0.6 
MeV. 

Our result is about 100% greater than Dashen's result. This result is in better 
agreement with earher references 0, |^, and with the recent investigation [Q, but 
in disagreement with Baur and Urech ||2|. Given the uncertainty of our result, we 
feel more comfortable by saying that we find a modification of Dashen's theorem of 
between 160% and 240% . 



6. Conclusions 



The calculation of nonleptonic amplitudes is in general one of the most difficult 
tasks for analytic strong interaction techniques. The elecromagnetic mass differences 
of the pseudoscalar mesons seems to us to be the most favorable case to attempt a 
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Table 2: Am€^' results. 



Narrow- width Pa{s)^ Pv i^) 







(MeV) 


(MeV) 


Subtr. 




2.56 


1.80 


Elastic 




0.92 


0.92 


Str. Fn. Ki int. 


St. 


0.05 


0.07 


Str. Fn. K* int. 


St. 


-0.18 


-0.18 


Total calculated 




3.35 


2.61 


Daslien 




1.27 


1.27 



controlled calculation. There turn out to be several favorable circumstances that help 
in this endeavor. As we have exploited above, the relevant current-current products 
have several connections to known phenomenology, and have important constraints 
due to the long distance chiral behavior and the short distance properties of QCD. 

The calculation of the known pion mass difference was quite successful. It turns 
out that intermediate mass scales (around 1 GeV) are the most important for this 
matrix element, and these are well represented by resonance contributions. In fact 
this structure is already visible in the old calculation in the soft pion limit given by 
Das et al. where the vector and axial vector spectral functions determine the mass 
difference in the chiral limit. There are calculable corrections and even new diagrams 
that come in as one includes a non-zero pion mass, but the pion mass is still small 
enough that one does not change the general anatomy of the matrix element. 

In the case of the kaon mass difference, the experimental result is not known. 
We find a large deviation from the prediction of Dashen's theorem, which is valid in 
the limit of massless kaons. While the magnitude of this effect is larger than most 
5'f/(3) breaking effects in chiral calculations, we stress that its origin is in reasonably 
well-known and mundane effects, and does not represent any breakdown of chiral 
symmetry. The main effect seems to be the kaon mass in the propagator of the Born 
diagram, which hence is a rather long distance effect, while the remaining dependence 
comes from the known shift in resonance masses due the the strange quark mass. This 
mass difference is important for the extraction of the u — d quark mass difference. 

Appendix 

The notation used for the 0{E'^) chiral terms in the Lagrangian of Eq. (|19D is 
U = exp (^iV2<^ / , 
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X 



V 

25o 



rriu 
md 
m, 



(122) 



where f^, a'^ are the external fields. In order to include EM one needs to define 



= , 



eQA 



(123) 



where is the photon field, and should not be confused with the axial-vector an- 
tisymmetric tensor field which has two Lorentz indices, A^^^. Q is the quark charge 
matrix, for the u, d and s quarks. 







Q 



2 
3 

-k 








-i 



(124) 



The notation used for Lagrangian containing the chiral couplings of the vector 
and axial-vector meson resonances, Eq. (^) is 



f: 



V. 



Ill/ 



A 



fll/ 



u 



R 



U 



v 



V2 

K*- 



"1 + A 



/ .0 

\/2 ' \/6 



V 



p 

V2 v/6 



*o 



^ + 4 ^? 



V2 V6 



*o 



K 



^1^ 



2 . 



(125) 



XJ^fiv 
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]-{u^[dx-iivx + ax)]u + u[dx-i{vx-ax)]u^} . (126) 



From the kinetic terms of the Lagrangian in equation (|T9D one derives the free 
propagator for the antisymmetric tensor field representation 



{Q\TR^,{x)Rp„{y)\Q) -- 



X 



M2 J 27r4 M^-k^- ie 
g^^pQua (m^ - A;^) + g^pkuk^ - g^ak^kp - {/i ^ u) 



(127) 



where the normahzation is given by 



{0\R^^\R{e,p)) = j^[pMp)-p,ep{p)] . (128) 
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